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При травмах нервно-мышечной системы конечностей для оценки тяжести патологического процесса и выявления динамики восстановления 
нервов в процессе курса лечения используют характеристики сигнала, регистрируемого при раздражении нерва электрическим током. В 
статье рассматривается его вейвлет-преобразование. Приводится анализ вейвлет-коэффициентов для различных групп пациентов. 


Введение 

Для исследования состояния нервно-мышечной 
системы используется метод электронейромиографии 
(ЭНМГ). ЭНМГ -метод основан на регистрации и анали- 
зе биоэлектрической активности мышечных и нервных 
волокон, как спонтанной, т.е. отражающей состояние их 
в покое и при произвольном мышечном напряжении, так 
и вызванном, т.е. обусловленной электрической стиму- 
ляцией нерва или мышцы различной интенсивности и 
частоты [1]. Электронейромиография позволяет судить 
об уровне поражения, о степени тяжести травмы. Она 
также используется в качестве оперативного контроля 
адекватности воздействия электростимуляции и для 
оценки скорости восстановительных процессов в не- 
рвно-мышечном аппарате. Современная аппаратура для 
проведения ЭНМГ позволяет регистрировать следую- 
щие ответы мышц: М-ответ, Н-рефлекс, Р-волну [1]. 
Среди ответов выделяют: 

- М-ответ - вызванный электрический ответ мышцы, 
возникающий при электрическом раздражении дви- 
гательных волокон нерва; 

- Н-рефлекс - рефлекторный, возникающий в мышце 
при электрическом раздражении чувствительных 
волокон нерва; 

- Р-волна, возникающая в мышце при электрической 
стимуляции двигательных аксонов нерва. 

М-ответ - это вызванный электрический ответ мыш- 
цы, являющийся суммарной одновременной реакцией 
двигательных единиц мышцы в ответ на электрическое 
раздражение нерва, график которого представлен на 
рис. 1 [1]. 

Форма, амплитуда и латентный период М-ответа оп- 
ределяют функциональное состояние нервных волокон 
и изменяются при травме. 

Построение математической модели, описывающей 
форму М-ответа, является актуальной задачей. Анализ 
параметров этой модели у различных групп больных 
позволяет произвести их группировку; оценить тя- 
жесть патологического процесса, дать оценку скорости 
восстановительных процессов в нервно-мышечном 
аппарате у индивида в курсе лечения после проведённой 
процедуры. 

В настоящее время для математического описания 
и анализа сигналов все больше используют вейвлеты. 
Термин “сигнал” обозначает набор численно зафикси- 
рованной информации о каком-либо процессе, объекте 
и т.п. Под “анализом” сигнала имеется в виду не только 


его математическое преобразование, но и получение 
на основе этого преобразования выводов о специфике 
соответствующего процесса или объекта. 

На основе понятия о векторном пространстве об- 
щепринятым подходом к анализу сигналов 8(4) стало 
их представление в виде взвешенной суммы простых 
составляющих - базисных Рзі-функций ѵ|/ к (4), помно- 
женных на коэффициенты С, [2]: 

Т)=Е Х - ■ч'-о. 


Так как базисные функции ѵ|/ к (1) представляют^ 
заданными как функции определенного вида, то только 
коэффициенты С к содержат информацию о конкретном 
сигнале. Таким образом, можно говорить о возможности 
представления произвольных сигналов на основе рядов 
(1) с различными базисными функциями. 

Вейвлет - это некоторая Рзі-функция, которая долж- 
на обладать следующими свойствами [2]: 

- График такой функции должен осциллировать вок- 
руг нуля в окрестности некоторой точки на оси I, 
причем -№ 

|ѵ( г ) 8т = о. 


- Норма этой функции должна быть конечной: 

/ -н® 

М|= МѴ(Т Т <+“■ 

Ѵ-со у 

Рзі-функция создается на основе той или иной ба- 
зисной функции \| / 0 (4), которая определяет тип вейвлета. 
Базисная функция должна удовлетворять всем требова- 
ниям, которые отмечены для Рзі-функции. Она должна 
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обеспечить выполнение двух операций: 

- смещение по оси времени 1 - \|/ 0 (1-Ь) при ЬеК; 

- масштабирование - а' 1/2 \|/ (1/а) при а>0 и аеК, где К 

- область ограничения сигнала 8(1:). 

Тогда Рзі-функцию можно записать: 

ч <( хл >“'' 2 


При заданных а и Ь функция \|/(1) и есть вейвлет. Вей- 
влеты являются вещественными функциями времени 
1, параметр Ь в (2) задает положение вейвлетов по оси 
времени, а параметр а - масштаб. (2) 

При вейвлет-преобразовании выбор типов базисных 
функций (вейвлетов) намного более обширен, чем при 
преобразовании Фурье. В качестве вейвлетов могут 
использоваться ортогональные и биортогональные непе- 
риодические функции, функции, имеющие глобальный 
экстремум и быстрое затухание на бесконечности и т.д. 
Все это дает обширные возможности для представления 
различных сигналов. 


Прямое непрерывное вейвлет-преобразование 
сигнала 8(1) задается путем вычисления вейвлет-коэф- 
фициентов по формуле [2]: 



где К - область ограничения сигнала 8(1), а,ЬеК, а^О, 
\|/(1) - вейвлет- функция, параметр Ь - задаёт положение 
вейвлет-функции, а параметр а - её масштаб. Вейвлет- 
коэффициенты определяются интегральным значением 
скалярного произведения сигнала на вейвлет-функцию 
заданного вида. 

Прямое вейвлет-преобразование можно рассматри- 
вать как разложение сигнала по всем возможным сдви- 
гам и растяжениям/сжатиям сигнала 8(1) или некоторой 
произвольной функции. При этом параметры а и Ь могут 
принимать любые значения в пределах указанных выше 
областей их определения. Заметим, что преобразование 
Фурье также можно рассматривать как разложения по 
сдвигам (имеется в виду фазовый сдвиг гармоник, зада- 



Рис. 2. Графики восстановленного при помощи вейвлет-преобразования сигнала М-ответа 
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Рис. 3. Области расположения вейвлет-коэффициентов 


ющих положение их графиков) и растяжением/сжатием 
(определяемым значением амплитуд гармоник), но 
применительно к одной функции (синусоиде), не очень 
удобной для представления локальных особенностей 
сигналов. 

Обратное непрерывное вейвлет-преобразование 
осуществляется по формуле [2]: 


Т)= 




5а 5р 
02 ? 


где К - константа, определяемая функцией \|/. 


Вейвлет-преобразование М-ответа 

Рассмотрим применение вейвлет-преобразо- 
вания для математического описания вызванного 
потенциала мышцы (М-ответа) с целью выявления 
его особенностей, которые характеризуют свойства 
нервно-мышечной системы человека. В частности, іШ 
терес представляет исследование соответствия между 
особенностями М-ответа и степенью тяжести травмы 
нерва конечности. 

Для исследования были взяты две группы людей: 
здоровых и больных, с различной степенью тяжести 
травмы. Для математического описания М-ответа в 
качестве базисной вейвлет-функции используется 



вейвлет “мексиканская шляпа”. Его временной образ 
описывается аналитическим выражением [3]: 

О т 2 \ е -і/2 
- )• ' ■ 

В соответствии с формулой (2) перейдем к вейв- 
лету: 


ч'СЛѴч'о 


/т (3 


1 - 


Р\ 2 ^ 


. ~ ^ К 


Прямое и обратное вейвлет-преобразование осу- 
ществляется по формулам (3) и (4). Все преобразова- 
ния проводились при помощи системы компьютерной 
математики МаіЬсасІ 2001. 

Для каждого индивида из представленных групп 
получены вейвлет-преобразование сигнала М-ответа 
и графическое изображение области расположения 
вейвлет-коэффициентов. Область расположения значи- 
мых (ненулевых) значений определяется параметрами 
Ь и а, где а - определяет масштаб амплитуды сигнала 
М-ответа, а Ь - расположение этой амплитуды по оси 
времени I. 

Анализ результатов приводит к следующим выво- 
дам: 

1. М-ответ у здоровых людей представляет собой 
сигнал, который имеет один положительный пик с 
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определенной амплитудой (рис. 2, а), а у больных 
- с несколькими пиками и амплитудами меньшими 
по величине (рис. 2, б, в). Это связано с тем, что 
сигнал М-ответа отражает одновременную реакцию 
двигательных единиц мышцы на электрическое раз- 
дражение нерва. У больных пациентов в результате 
травмы двигательные единицы мышцы реагируют 
неодновременно, и это отражается в виде нескольких 
пиков. 

2. У здоровых людей (рис. 3, а) области локали- 
зации значимых вейвлет-коэффициентов отли- 
чаются от областей больных пациентов (рис. 3, 
б, в) как характером расположения (по параметру 
а), так и их числом (по параметру Ь). Количество 
положительных пикой сигнала М-ответа связано 
с количеством областей локализации значимых 
вейвлет-коэффициентов. 

3 . В группе больных пациентов по числу расположения 
вейвлет-коэффициентов (по параметру Ь) выдели- 
лось две подгруппы. Анализ состояния пациентов, 
проведенный врачом-экспертом, позволяет сделать 
вывод, что у пациентов с легкой степенью тяжести 
травмы области локализации две (рис. 3, б), а с 
тяжелой - три (рис. 3, в). 

4. Границы (по параметру а) локализации значимых 
вейвлет-коэффициентов для различных групп 
пациентов также отличаются. Так у здоровых они 


расположены в пределах от 5 до 25, у больных с 
легкой степенью тяжести травмы - от 4 до 17, с 
тяжелой - от 3 до 10. 

Таким образом, использование вейвлет-преобразова- 
ния для математического описания М-ответа позволяет 
восстановить сигнал, а по характеру и количеству об- 
ластей расположения значимых вейвлет-коэффициентов 
можно оценить степень тяжести травмы. 

В дальнейшем, используя математическое описание 
М-ответа, предполагается построить обобщенные по- 
казатели, с помощью которых можно было бы оценить 
тяжесть патологического процесса и скорость восста- 
новительных процессов в нервно-мышечном аппарате 
в курсе лечения. 
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